 figure(11);
 count = 0; 
 for h=[1 4 8 12]
    count = count+1;
 subplot(2,2,count);

plot(nn(floor(win/2)+h:floor(win/2)+size(MSFE_var{h}+h,2)+h-1),[ones(size(MSFE_var{h}+h,2),1) prctile(MSFE_var{h}./MSFE_ar{h},[16 50 84],1)'],'r');datetick('x','keepticks')
set(gca, 'XTick', nn(17:80:end));
 anni=datevec(nn(16+1:80:end,:));
set(gca, 'XTickLabel', anni(:,1));  

hold on
plot(nn(floor(win/2)+h:floor(win/2)+size(MSFE_var{h}+h,2)+h-1),[ones(size(MSFE_var{h}+h,2),1) prctile(MSFE_var{h}./MSFE_ar{h},[50 ],1)'],'k'); datetick('x','keepticks')
hold on
plot([datenum(1914,9,15) datenum(1914,9,15)], [0 9]); axis([696410 732021 0 9] )
plot([datenum(1931,6,15) datenum(1931,6,15)], [0 9]); axis([696410 732021 0 9] )
plot([datenum(1941,9,15) datenum(1941,9,15)], [0 9]); axis([696410 732021 0 9] )
plot([datenum(1951,3,15) datenum(1951,3,15)], [0 9]); axis([696410 732021 0 9] )
plot([datenum(1971,9,15) datenum(1971,9,15)], [0 9]); axis([696410 732021 0 9] )
plot([datenum(1983,12,15) datenum(1983,12,15)], [0 9]); axis([696410 732021 0 9] )
hold off
title('12-QUARTERS HORIZON')
% set(gcf,'DefaultAxesColorOrder',[0 0 0;1 0 0;1 0 0;0 0 1]); 
end
 